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Under  the  sponsorship  of  the  ONR  Contract  N00014-86-K-0533  we  have  published 
3  referenced  journal  and  conference  papers  and  5  technical  reports  on  time  domain  elec¬ 
tromagnetic  waves  in  multilayer  media. 
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'  The  analysis  of  resonance,  input  impedance  and  radiation  of  the  elliptic  disk,  mi¬ 
crostrip  structure  is  rigorously  formulated  in  this  paper,  using  the  Scalar  and  Vector 
Mathieu  Transforms.  With  the  help  of  these  transforms,  the  resonance  frequencies  of 
the  structure  can  be  derived  exactly  using  Galerkin’s  method  and  approximately  using  a 
perturbations!  approach.  Expressions  for  the  input  impedance  and  the  radiation  pattern 
are  also  obtained. 

Theory  for  quasi-TEM  modes  propagating  in  a  transversely  inhomogeneous  (multi¬ 
dielectric)  longitudinally  uniform  transmission  line,  previously  derived  for  time-harmonic 
waves,  is  derived  for  transient  signals.  It  is  seen  that,  while  the  starting  point  for  the  theory 
is  completely  different,  the  result  is  similar  to  the  time- harmonic  theory,  and  previously 
derived  properties  for  propagating  modes  also  apply  in  the  transient  case.  The  range  of 
applicability  is  discussed  with  a  simple  example. 

Exact  image  method,  recently  introduced  for  the  solution  of  electromagnetic  field 
problems  involving  sources  above  a  planar  interface  between  two  homogeneous  media,  is 
shown  to  be  valid  also  for  sources  located  in  complex  space,  which  makes  its  application 


possible  for  Gaussian  beam  analysis.  It  is  demonstrated  that  the  Goos-Hanchen  shift 
and  the  angular  shift  of  a  TE  polarized  beam  are  correctly  given  as  asymptotic  results 
by  the  exact  reflection  image  theory.  Also,  the  apparent  image  location  giving  the  correct 
Gaussian  beam  transmitted  through  the  interface  is  obtained  as  another  asymptotic  check. 
The  present  theory  makes  it  possible  to  calculate  the  exact  coupling  from  the  Gaussian 
beam  to  the  reflected  and  refracted  beams  as  well  as  to  the  surface  wave. 

Exact  image  method,  recently  introduced  for  the  solution  of  electromagnetic  field 
problems  involving  sources  above  a  planar  interface  between  two  homogeneous  media,  was 
originally  applied  to  a  restricted  pair  of  medium  parameters  to  obtain  a  well  behaved  image 
located  in  the  ’proper’  complex  half  space.  It  is  demonstrated  here  with  an  example  that 
the  ’proper’  half  space  limitation  can  be  released  to  increase  the  applicability  of  the  exact 
image  theory.  However,  it  is  shown  that  for  certain  media,  numerical  difficulties  in  the 
field  integrals  may  be  encountered,  due  to  crossing  of  branch  cut  lines  on  the  complex 
integration  plane.  This  may  occur  when  the  medium  where  the  field  is  to  be  calculated 
is  more  lossy  than  the  other  medium.  Methods  for  the  image  integration  to  obtain  better 
convergence  for  more  general  media,  are  discussed. 

Simple  approximation  for  diffraction  surface  currents  on  a  conducting  half  plane,  due 
to  an  incoming  plane  wave,  is  given  in  terms  of  a  line  current  (monofile)  in  complex  space. 
When  compared  to  the  approximation  by  a  current  located  at  the  edge,  the  diffraction 
pattern  is  seen  to  improve  by  an  order  of  magnitude  for  a  minimal  increase  in  computation 
effort.  Thus,  the  inconvenient  Fresnel  integral  functions  can  be  avoided  in  quick  calcula¬ 
tions  of  diffracted  fields  and  the  accuracy  is  seen  to  be  good  in  other  directions  than  along 
the  half  plane.  The  method  can  be  generally  applied  to  problems  involving  planar  metal 
edges. 

Many  integrated  circuits  contain  strip  lines  at  different  heights  that  run  parallel  or 
perpendicular  to  each  other.  And  we  have  investigated  reliable  models  for  these  struc- 


tures.  First  the  capacitances  associated  with  two  offset  parallel  strips  at  different  heights 
between  ground  planes  are  computed  using  the  conformal  mapping  approach.  As  an  ex¬ 
tension,  a  simplified  circuit  of  parallel- plate  lines  with  transverse  ridges  is  introduced  to 
model  two  parallel  strips  with  perpendicularly  crossing  strips  on  top.  We  treated  it  as 
a  distributed  circuit  consisting  of  transmission  lines  segments  with  periodical  capacitive 
loading.  In  order  to  calculate  the  coupling  between  two  lines,  we  reduced  this  structure 
to  two  equivalent  single  line  circuits,  viz.  the  even  mode  and  the  odd  mode  circuits.  The 
Laplace  transform  approach  can  be  easily  applied  to  find  out  the  transient  response.  The 
numerical  computation  carried  out  for  various  environments  shows  that  the  crossing  strips 
will  cause  serious  trouble  for  signals  with  a  rise  time  of  less  than  50ps  to  propagate  along 
a  distances  of  2cm  or  longer. 

We  have  also  proposed  a  way  to  account  for  the  effect  of  complicated  geometry  from 
the  point  of  view  of  continuous  transmission  line  model  by  letting  the  coupling  between 
parallel  lines  in  multilayered  integrated  circuits  to  be  nonuniform.  Previously,  we  devised 
a  scheme  that  combines  the  method  of  characteristics  and  perturbational  series  to  simplify 
the  computation  of  the  transient  response  from  the  coupled  transmission  line  equations.  A 
new  transformation  for  decoupling  now  enables  us  to  generalize  this  formulation  to  calcu¬ 
late  the  near-end  and  far-end  crosstalks  to  very  high  accuracy,  given  arbitrary  positional 
dependence  for  both  capacitive  and  inductive  coupling  coefficients. 

The  transient  response  of  fundamental  sources,  such  as  dipole  and  line  current,  was 
carefully  analyzed.  With  the  double-deformation  technique,  which  is  a  modal  technique 
based  on  identification  of  singularities  in  the  complex  frequency  and  wave  number  planes, 
we  are  able  to  obtain  both  early  and  late  time  response  very  efficiently.  Some  results  for 
vertical  electric  dipole  excitation  on  a  two-layer  medium  have  been  published.  Recently, 
we  have  discovered  a  general  scheme  of  breaking  up  the  integrands  so  that  sources  with 


arbitrary  time  and  space  dependence  can  be  incorporated  into  our  formulation  without 
sacrificing  convergence. 

Basically,  vias  in  a  multilayered  integrated  circuits  are  treated  like  transmission 
lines  with  loadings  where  they  encounter  holes  in  ground  planes  separating  different  lay¬ 
ers.  We  have  modeled  a  ground  plane  with  a  hole  and  a  circular  conductor  at  the  center 
of  the  hole  as  a  radial  waveguide,  which  in  turn  is  connected  to  the  via,  another  section 
of  transmission  line.  Thus  by  computing  the  characteristic  impedance  of  the  former,  we 
have  derived  the  equivalent  load  impedance  of  the  via  hole.  The  load  impedance  is  one 
important  parameter  in  calculating  the  transient  propagation  along  vias. 
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